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Эволюционное моделирование задач 
многокритериальной оптимизации в системе 
исследовательского проектирования 
морских судов 


Рассматривается подход и концептуальный базис эволюционного моделирования на современном этапе 
развития учения о естественной генетике в задачах исследовательского проектирования. Сформулированы 
теоретические принципы формализации проблемы с использованием достижений интеллектуальных тех- 
нологий и высокопроизводительных средств вычислений. Приведены примеры реализации разработанного 
подхода в практических приложениях исследовательского проектирования морских судов. 


Парадигма эволюционного моделирования в системе 
исследовательского проектирования морских судов 


Создание формальной системы эволюционного моделирования связано с решением 
задач многокритериальной оптимизации, описывающих поведение судна в сложной 
динамической среде с помощью нового подхода к интерпретации динамики внешней 
среды, осуществляется на основе реальных данных о климатических спектрах морского 
волнения [1]. В отличие от имеющихся исследований [2-8] разработанная концептуаль- 
ная модель эволюционного моделирования описывает существенные свойства ис- 
следуемой проблемной области в рамках нового подхода к решению задач многокри- 
териальной оптимизации [9-12]. Функциональный блок эаволюционного моделирования 
определяет построение и интерпретацию математической модели на базе интеллектуаль- 
ной системы (ИС) исследовательского проектирования. Решение задач многокритериаль- 
ной оптимизации ведется в рамках формального и строгого определения понятий и 
закономерностей естественной генетики [13] и формального аппарата теории систем [14]. 

Проблема оптимизации в задачах моделирования. Проблема оптимизации при мо- 
делировании поведения нелинейных динамических систем [15-17], особенно неста- 
ционарных систем, функционирующих в сложной динамической среде [15], [18], пред- 
ставляет большой теоретический и практический интерес в связи с созданием новой 
техники и технологий на базе достижений современной компьютерной математики и 
методов искусственного интеллекта [19]. Среди исследований в области оптимизации 
сложных динамических систем на современном этапе можно выделить следующие 
направления исследований [9]: 

— использование классических методов и моделей, адаптированных с учетом реаль- 
ных процессов взаимодействия объектов моделирования с внешней средой при различном 
уровне возмущающих воздействий; 

— использование нетрадиционных подходов к оптимизации, основанных на дости- 
жениях интеллектуальных технологий и высокопроизводительных средств вычислений. 
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Оба эти направления связаны с решением задач оптимизации в сложной динами- 
ческой среде, которая характеризуется НЕ-факторами [20]: неопределенность, неполнота, 
нелинейность, нестационарность. В этих условиях нельзя сразу отдать предпочтение 
какой-либо модели и использовать принцип конкуренции [18], позволяющий на основе 
генерации альтернатив выбирать предпочтительное решение. Принцип конкуренции 
сформулирован в рамках концепции мягких вычислений [21], объединяющей нечеткий 
логический базис с нейросетевыми моделями и генетическим алгоритмом (рис. 1). 

Реализация принципа конкуренции [18] предусматривает использование критериаль- 
ной базы оценки вариантов решения текущей задачи моделирования, сгенерированных в 
процессе эволюции объекта моделирования. Среди методов и моделей, используемых 
при реализации концепции мягких вычислений, следует выделить подход Беллмана — 
Заде [22], метод иерархий Саати [23], а также более простые оценки, основанные на 
модели выбора [24], функции выбора [25] и методе эталонов [26]. 

Поток информации, реализующий концептуальную модель эволюционного моде- 
лирования в системе исследовательского проектирования при интерпретации решений 
на основе интеллектуальных технологий и высокопроизводительных вычислительных 
средств, представлен на рис. 1. 


Поток информации в мультипроцессорной вычислительной среде 


Реализация конкурирующей 
Принцип вычислительной технологии: 


конкуренции параллельная обработка информации на базе 
традиционных и нейросетевых моделей 


Принцип Построение адаптивных моделей знаний и 
формализации Г реализация механизма логического вывода в 
нечетких знаний мультипроцессорной вычислительной среде 


Рисунок 1 — Организация параллельной обработки информации 
при функционировании системы исследовательского проектирования 
морских судов 


Исследования в области оптимизации показывают, что достаточно эффективным 
направлением, реализующим новые подходы к этой проблеме, являются эволюцион- 
ные модели, основанные на использовании генетического алгоритма [10-12], [19]. Суть 
эволюционного подхода сводится к следующему: 

— Фиксируется множество объектов Х, параметры которых связаны друг с другом 
посредством структуры, и необходимо выбрать наилучшие структуры на основе за- 
данного критерия оптимальности, который формируется на основе свойств объектов и 
может быть представлен в виде различных моделей, причем каждому объекту мно- 
жества Х сопоставляется определенное значение критерия Е(Х). 

— Множество объектов произвольно, и необходимо построить представление $ ис- 
ходного множества объектов в другом множестве, обладающем некоторой структурой, 
в том числе и векторным пространством, причем связь между исследуемыми 
объектами описывается представлением ф: Х->5, а объекты выступают в качестве 
потенциальных решений задачи поиска экстремума с помощью поискового алгоритма. 
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— Процесс оптимизации на основе эволюционного моделирования состоит в по- 
строении множества объектов-решений, для которых выполняется следующее условие: 


Х* = аго тах Е [ф \($]]. (1) 
ЗЕЗ 

Алгоритмы оптимизации представляют собой дискретные процедуры, состоящие в 
преобразовании исходного множества решений путем отбора наилучших по критерию Е 
(функции пригодности) решений и генетических операторов, основными из которых 
являются кроссовер и мутация. 

Функциональные блоки, с помощью которых реализуется процесс обработки ин- 
формации в задачах эволюционного моделирования, представлены на схеме рис. 2. 
Здесь выделена последовательность преобразования информации многокритериальной 
оптимизации — от базы данных ИС исследовательского проектирования — до выдачи 
практических рекомендаций для принятия обоснованного решения. 


Информационная 
среда анализа и 
интерпретации 

принятия решений 


Компоненты 
информационной 
поддержки 


Функциональные блоки о =. 
исследовательского данных 


проектирования 


Рисунок 2 — Функциональные блоки системы интеллектуальной поддержки 
конструктора, реализующие процедуру эволюционного моделирования 
на базе технологии интеллектуальных систем 


Ниже дается краткая характеристика принципов, при построении функциональ- 
ного блока эволюционного моделирования, использование которых ориентировано на 
новые поколения систем исследовательского проектирования морских судов с учетом 


реальных данных о климатических спектрах морского волнения в районе их эксплуа- 
тации [15], [27]. 


Принципы обработки информации, 
обеспечивающие быстродействие алгоритмов 
эволюционного моделирования в задачах 
исследовательского проектирования морских судов 


Повышение эффективности функционирования, достоверности оценки и прогноза 
исследуемой ситуации в вычислительном комплексе эволюционного моделировния, 
использующих интеллектуальные технологии, достигается в рамках парадигмы обра- 
ботки информации в мультипроцессорной вычислительной среде (рис. 3) [18]. Расши- 
рение функциональных возможностей системы обеспечивается на основе принципов 
открытости, сложности и конкуренции в рамках концепции мягких вычислений [15], [21]. 
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состояния объекта объект 


ивнешней среды 
ве денных данных и _ 
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< онетииюинныи принадлежности 


Рисунок 3 — Поток ом при реализации концепции мягких вычислений 


Реализация принципа конкуренции [18] предусматривает использование кри- 
териальной базы оценки вариантов решения текущей задачи моделирования, сгенери- 
рованных в процессе эволюции объекта моделирования. Среди методов и моделей, ис- 
пользуемых при реализации концепции мягких вычислений, следует выделить подход 
Беллмана — Заде [22], метод иерархий Саати [23]. Необходимо отметить, что подход 
Беллмана — Заде и метод анализа иерархий реализуют алгоритм многокритериальной 
оптимизации в нечеткой среде. Причем первый из них основан на построении экспертной 
матрицы критериев и ограничений с последующей операцией минимаксной композиции, 
тогда как метод анализа иерархий предполагает построение иерархической структуры в 
виде матрицы парных сравнений. 

Решение задачи многокритериальной оптимизации в рамках подхода Беллмана — 
Заде состоит в достижении цели при заданных нечетких ограничениях. Пусть некото- 
рая альтернатива х обеспечивает достижение цели со степенью ис(х) и удовлетворяет 
ограничениям со степенью ис(х). Тогда степень принадлежности альтернативы х ре- 
шению задачи равна минимальному из этих чисел, т.е. (функция принадлежности) мо- 
жет быть интерпретирована в зависимости от числа целей 7; и ограничений у; целей, а 
также их важности: 


но(х) = пи {ис(х), ис(х}. 
ис(х) = пи {с1(Х), ... ‚ Нат(Х), иси(Х), ... › Ис}. (3) 
ио(х) = тш {Аиси(Х), ... › АИаь(Х), УшскХ), ... ‚ УтИст(Х)}. 
С учетом введенных обозначений функция ир(х) выражается как 


1 м. Ис(У) при хЕ р 


(4) 
1прихе Х/Х°. 


Ир(х) = 


В соответствии с определением нечеткое решение рассматриваемой задачи имеет 
ВИД: 
ыо(х) = пи {ыс(х), ис()}; 
тие), Е, Ис(у) } при хе Х°; 


Ир(х) = м (5) 
ис(х) прихе Х/хХ°. 


Метод анализа иерархий является одним из наиболее распространенных при ин- 
терпретации данных в ИС [23]. В отличие от подхода Беллмана — Заде, метод иерархий 
предполагает представление исходной информации в виде относительных весов приз- 
наков, определяющих текущую ситуацию. Суть этого метода состоит в реализации 
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последовательности операций анализа иерархической структуры. Выделенные в этой 
структуре элементарные подзадачи анализируются методом иерархий [23], с помощью 
которого выделяются приоритеты (веса) признаков в целях выбора наилучших из 
них с использованием алгебраической теории матриц и экспертных процедур: 


И =А ал; 
Пи * Лш 81 (6) 
Е Е В м. 
Пи р Л пт 8т 


где \\ — обратно симметричная матрица значений парных сравнений признаков от- 
носительно данного атрибута; л-— нормированный вектор весов признаков; Аах — наи- 
большее собственное значение матрицы \№; П - результат определения глобальных 
приоритетов признаков Ш, ..., Пм; М — число признаков; ль ( 11....п, ]=1,...ла) — 
относительный вес 1-го признака по }-му атрибуту; 2; — относительный вес ]-го атрибута. 

Матрица (6) обладает свойством 

П-М=0, (7) 
где [ — единичная матрица. 

Контроль надежности выводов, полученных при использовании метода анализа 
иерархий в задачах принятия решений реализуется с использованием индекса согласован- 
ности /[х. Величина [5 определяется для каждой матрицы парных сравнений по вектору 
приоритетов Г. Вектор Г вычисляется как собственный вектор матрицы парных 
сравнений. 

Процедура оценки индекса согласованности [5 осуществляется по соответст- 
вующему максимальному собственному числу АД: 

15 =(#-п)/п-1, (8) 

где п > 1 — число сравниваемых элементов (признаков или альтернатив), т.е. и> 1. 
Для обратносимметричных матриц, к которым относятся матрица парных сравнений (6), 
й> п ‚ поэтому индекс согласованности — величина неотрицательная и конечная. 

Полученные индексы согласованности сравниваются с соответствующими сред- 


ними Вз= В (п), которые определяются при случайном выборе решения. Эти значения 
приведены в табл. 1. 


Таблица 1 — Критические значения индекса согласованности 
п 2 3 4 5 6 7 8 9 
В$ (и) 0,1 10,58 | 0,90 | 1,12 | 1,24 | 1,32 | 1,41 11,45 


Оценкой надежности метода парных сравнений является отношение согласован- 


ности Ду: 
45 == [5 / Кс . (9) 
Если величина 45 не превышает 0,2, то полученные результаты являются до- 
стоверными. 


Особенность этого метода состоит в том, что результаты парных сравнений для 
двух сложных объектов достаточно сложно передать в виде точных цифр. Назначение 
цифр нередко кажется искусственным и достаточно произвольным. Действительно, 
представление результата сравнения в виде разумных цифр требует глубокого по- 
нимания задачи — насколько свойства рассматриваемых объектов влияют на дости- 
жение поставленной цели. 
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При анализе альтернатив в сложных ситуациях находят применение нечеткие 
оценки и граф-интерпретации, позволяющие рассматривать формирование действий 
конструктора в виде комбинаторных задач на графах. Область допустимых решений 
в этом случае определяется уровнем нечеткости графа, соответствующим нечеткой 
оценке. Задаваясь такой оценкой, можно определить область допустимых решений, 
а затем в этой области вести поиск оптимальных решений на графах [15]. 

Путь на графе для конкретной ситуации объединяет совокупность вершин с опре- 
деленным порядковым номером. С каждым таким путем сопоставляются оценки, опреде- 
ляющие риск и временные затраты. Обозначив $ = (51,...,3т) — множество всех путей в 
граф-модели плана реализации решения, можно установить экстремальный путь, опреде- 
ляемый соотношением: 


5*= шаех м шш{нио($1),..., Ис(Зв), (10) 
где ис($1),..., ис(Зв) — функция принадлежности множества путей. 


Концепция эволюционного моделирования в системе 
исследовательского проектирования 


Формальное описание общей модели эволюционного моделирования реализуется с 
помощью генетического алгоритма СА [11], [12]: 

СА = (Ро, Г, 15, В, Е, 0, (11) 

где Ро = (рь,...,рт,); р. — решение задачи, представленное в виде хромосомы, кодиро- 

вание которой осуществляется с помощью кода Грея [16]; Г. — целое число (размер 

популяции); [| — целое число (длина каждой хромосомы в популяции); В — отображение, 

определяющее рекомбинацию (кроссовер (сгоззоуег), мутация (ти@йоп); Е — функция 
оптимальности (В пе5$ Кисйоп); {— критерий остановки. 

При использовании СА в задачах оптимизации простой рандомизированный 
перебор при поиске минимума по методу Монте-Карло можно заменить на хранение 
популяции лучших с точки зрения функционала объектов и добавления новых объектов 
посредством мутации одного объекта, либо посредством кроссовера — обмена подоб- 
ными частями у двух случайных объектов. Мутация — это почти полный аналог случай- 
ного шага в методе Монте-Карло, а кроссовер и популяция привлечены из биологии, 
где естественный отбор оптимизирует генотип [13]. Поток информации при реали- 
зации СА в задаче оптимизации представлен на рис. 4. 


Ро и(х)|1....0)}. (12) 


Формирование исходной базы знаний 
(популяции) из т-структур 


Вычисление меры полезности 
для каждой структуры 


Вычисление вероятности выбора 
каждой структуры 


Формирование новой базы знаний на основании 
вероятности выбора и применения генетических 
операторов к отобранным структурам 


Рисунок 4 — Парадигма преобразования информации при реализации СА 
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Современные исследования СА в основном связаны с моделированием как средст- 
вом изучения поведения алгоритмов и придания им новых свойств. Структурно-ин- 
вариантный анализ позволяет применить к моделям СА практически тот же метод 
декомпозиции, что и для динамических систем. 

Из рисунка видно, что поиск более совершенных структур достигается с помощью 
СА, преобразующего текущую популяцию (Б3). Это преобразование обеспечивается за 
счет селективного смещения в пользу более устойчивых структур. Таким образом, цель 
СА состоит в абстрактном и формальном объяснении процессов эволюции и адаптации, 
а также создания аппарата проектирования искусственных программных систем, 
содержащих механизмы эволюции естественных систем. 

СА отличаются от других оптимизационных и поисковых процедур. Это от- 
личие состоит в следующем [11], [12], [19]: 
— работа не с параметрами, а с закодированными множествами; 
— осуществление поиска из популяции точек, а не из единственной точки; 
— использование для оценки информации целевой функции, а не ее приращений; 
— использование вероятностных, а не детерминированных правил. 

Основными генетическими операторами в задачах эволюционного моделирова- 
ния являются кроссовер и мутация. Кроссовер (скрещивание) моделирует передачу 
наследственности хромосомами и обуславливает целенаправленное приближение хро- 
мосом к оптимальному решению. Однако использование классических операторов скре- 
щивания (одноточечных или многоточечных кроссоверов) в отдельных случаях является 
малоэффективным. В работе [10] предложены другие операторы рекомбинации (рис. 5), 
уменьшающие число переборов и время сходимости алгоритма оптимизации. 


Рисунок 5 — Оператор рекомбинации: а: ; с; 4; — значения генов 


Одна из проблем эволюционного моделирования состоит в построении решения, 
позволяющего избежать попадания в локальный экстремум. Эффективная процедура 
выхода из локального экстремума связана с хромосомной мутацией (хромосомной 
перестройкой) [10]: дефишенси (концевые нехватки хромосом), дупликация (удвоение 
участков хромосом), инверсия (перестройка части генов в обратном порядке), транс- 
локация — перенос части генов в той же хромосоме на новое место (образование изо- 
хромосомы) (рис. 6). 
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Транслокация Результат 
хромосомы транслокации 


Рисунок 6 — Оператор транслокации 


Данные типы перестроек играют важную роль в естественной эволюции, поз- 
воляя переходить к другим видам [10]. Построенные на их основе операторы мутации 
СА позволяет обходить локальные экстремумы. 

Комбинация операторов кроссовера и мутации образуют схему замещения М в 
виде вектор-функции М(Р), компонентами которой являются квадратичные формы: 

М(Р) =<р°"Мр”, [5(ОРОГМ 5(2)р°....., 
[59р°ТМ $9"... ТМ $(1)р'>, (13) 
где М - матрица замещения; 5(К) — матрица перестановок, зависящие от элемента 
популяции К. 

Групповая структура на множестве Х коммутирует со схемой замещения и 
задается матрицами 5(К) в виде линейной группы матриц. Инвариантность схемы 
замещения относительно группы определяет симметрию СА. 

Другим эффективным направлением в совершенствовании методов оптимизации 
сложных систем, функционирующих в нестационарной динамической среде, является 
использование преимуществ нейросетевого логического базиса. Реализация этих пре- 
имуществ привела к разработке теории нейро-эволюционного моделирования, соче- 
тающего массовый параллелизм и нелинейное преобразование информации нейронных 
сетей с вычислительной технологией эволюционного моделирования. В рамках такой 
интеграции находят применение ансамбли нейронных сетей, а также интегрированные 
нейронные сети, использующие преимущества КВЕ-сетей, многослойного персептрона 
и самоорганизующихся карт Кохонена [28]. 


Последовательность операций, обеспечивающих 
процедуры эволюционного моделирования 


Рассмотрим особенности задачи эволюционного моделирования. Пусть Х — 
некоторое множество объектов (особей), среди которых надо найти оптимальный 
объект, а $ — множество двоичных строк составленных из 0 или 1. Заданы отобра- 
жения 

Е[Х К, Е2:Х55 и ЕЗ:$=>Х, (14) 
где 5 — множество действительных чисел. Функцию Е1 в СА называют функцией 
приспособленности. Задан оператор 

М1:$ 55, (15) 


переводящий $1 е Зв $2е 5 путем инверсии в выбранных позициях. 
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Выбор позиций, подвергаемых инверсии, недетерминирован. 
Введем оператор мутации 


мМ2:Х5Х, (16) 
определяемый как 
М2()ЕЕЗ(М КЕ2(Х))), хЕХ. (17) 


Операцию кроссовера определим следующим образом. Пусть $1ЕЗ в $2еЗ, аР 
множество номеров позиций, в которых строки различаются. Порождаемая при вы- 
полнении кроссовера строка обладает следующими свойствами: 

— строка может совпадать со строками $1 и $2 в позициях, не вошедших вр; 

— приблизительно половина позиций, вошедших в р, совпадает с позициями в 

строке $1, а другая половина позиций из Р совпадает с позициями в строке $2. 
Требуется найти хе Х такой, что 


Е (хутах {Е 1(7)| 2Х}, (18) 


т.е. ставится задача поиска строки с максимальным значением функции приспособлен- 
НОСТИ. 


Применяя последовательность «отбор — мутация» (рис. 7), можно направлять 
эволюцию элементов популяции к наиболее хорошим точкам пространства поиска. 


Формирование Отбор хромосом для 
начальной популяции репликации 


У 


Вычисление мер Операция 
пригодности пересечения 
Определение элитной 
Мутгация 


хромюсомы 


У 


Условие окончания об Е 
поиска аи 
Вывод, Формирование новой 
решения популяции 


Рисунок 7 — Обобщенная блок-схема генетического алгоритма 


При этом отбор производит селекцию подходящих элементов, а кроссовер и му- 
тация накладывают стохастический шум на процесс эволюции, осуществляя переход к 
исследованию новых областей и точек пространства. 


Нечеткая формальная система, обеспечивающая 
процедуры эволюционного моделирования 


Композиционная модель нечеткого вывода в динамической базе знаний ИС ис- 
следовательского проектирования описывает связь всех возможных состояний логи- 
ческой системы с управляющими воздействиями и формально записывается в виде 
(Х,В,У), где Х={хь,..., хи}; У{уь,....Ун! — базовые множества, на которых заданы входы 
А; и выходы В; нечеткой системы; К - нечеткое соответствие «вход — выход», которое 
строится на основе формализации нечетких стратегий. Принятие решений при функ- 
ционировании нечеткой системы знаний осуществляется на основе модели [15]: 


«ША Феп Ву, е1зе МА› Феп В....., е15е Е Ах Феп Вх». (19) 
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Здесь Ат,...,Ам — нечеткие подмножества, определенные на базовом множестве Х, 
а В.....,Вм — нечеткие подмножества из базового множества У. Эти нечеткие подмно- 


жества задаются с помощью функций принадлежности ид(х), ивкУ). (хех, у&). 


В=Ахв= {иль} (20) 


1=1 75 
ИЛИ 


Ив(х, у) = ид(х) ^ Ив(у) = пиш(иа(х), Ш (У)), (21) 
где нА(х), ив(у) — ФП элементов х, у множествам А и В. 
В качестве нечеткой импликации при реализации механизма логического вы- 
вода в ИС использованы алгоритмы Мамдани и Сугено [15], [29]. 
Алгоритмы преобразования информации при реализации механизма логического 
вывода в ИС представлены на основе продукционных моделей [15], [27]: 


< Х, еп У, е[зе 7>; (22) 
(Х >У(2)), 
где Х, У, / — нечеткие множества, определенные на универсальных мно- 


жествах Ч, У, \; множества Х, У, 7, интерпретируются как нечеткий вход и нечет- 
кий выход некоторой системы, отношения между которыми определяются моделью 
нечеткой системы. 

При интерпретации модели (22) используются различные встроенные процедуры, 
основанные на применении аппарата статистического анализа, методов математического 
моделирования, методов качественного исследования дифференциальных уравнений на 
фазовой плоскости и др. Наряду с методами традиционной математики и аппаратом не- 
четкой логики здесь находят также широкое применение теория ИНС [28] и СА [11], [12]. 

Механизм, обеспечивающий логику функционирования и выработку решений 
по мере поступления информации о параметрах ДО и внешней среды, основан на 
выводе решений с помощью формальной процедуры [15], [27]: 

Е; : Эк (В) > Ц; (23) 
(К =1,...0;]=1,....51=1,...Ю), 

где 5к(Е) — ситуация в момент времени & ; Ч;— решение; Е; — множество правил. 

Практическое приложение. Рассмотрим применение нейро-эволюционного мо- 
делирования в задаче контроля динамики сложного объекта с использованием концепции 
нейро-нечеткой системы АМЕБ [30]. Потенциально настраиваемыми параметрами моде- 
ли нечеткого вывода по Сугено в базе знаний адаптивной ИС, реализующей поддержку 
нейро-эволюционного моделирования, являются параметры функций принадлежности и 
коэффициенты в функциях, вычисляющих выход для нечетких логических правил. 

Необходимым условием применения СА является представление параметров за- 
дачи в форме хромосом, т.е. в виде наборов ген, которые могут принимать различные 
значения (аллели). В данном случае параметрами задачи является набор коэффициен- 
тов системы нечеткого вывода: р, 4; и г; [16]. Для характеристики параметров системы 
применялось двоичное кодирование. Каждый параметр кодировался в виде 16 двоич- 
ных значений. Определялась допустимая область значений каждого коэффициента. Затем 
данная область разбивалась на 2'° равных интервала, каждому из которых присваивался 
номер для кодирования в виде генов с применением кодов Грея. Длина получившейся 
хромосомы равна 16хМ, где М — число кодируемых параметров системы. 

В задаче подбора параметров нечеткой системы шаги СА несколько отличались 
от шагов классической реализации. Рассмотрим эти шаги. Инициализация заключается 
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в случайном выборе заданного количества хромосом. Исходная популяция была по- 
лучена посредством случайной модификации параметров системы, определенной ме- 
тодом экспертной оценки и кодирования их в виде хромосом. Начальная популяция, 
так же как и все последующие, представлена 100 хромосомами. Селекция хромосом 
заключалась в выборе тех хромосом, которые будут участвовать в создании потомков 
для следующей популяции. 

При разработке системы управления сложным ДО в нестационарной динами- 
ческой среде построение такой функции является нетривиальной задачей. Отбор хро- 
мосом осуществлялся по следующему принципу. Выбирается случайная хромосома и 
преобразуется в набор параметров системы нечеткого вывода. После чего производится 
имитационное моделирование поведения ДО. При этом в качестве параметров системы 
управления используется набор параметров, полученных из тестируемой хромосомы. 
Если параметры отвечают заданным критериям — хромосома отбирается для создания 
следующей популяции, в противном случае тестированию подвергается следующая слу- 
чайная хромосома. Такая селекция хромосом полностью отвечает принципу естествен- 
ного отбора и легко реализуема при наличии модели процесса. Для получения более 
адаптированных к различным условиям популяции хромосом имитационное моделиро- 
вание каждый раз проводилось с различными внешними условиями и характеристи- 
ками ДО. С целью уменьшения вероятности случайного отбора «плохой» хромосомы 
имитационное моделирование можно проводить несколько раз для одной хромосомы. 

На рис. 8 приведен график, отображающий результат уменьшения количества 


неудачных решений для каждой следующей популяции. 
80 


1 $ 9 1317212529 33374145 49 535761 65 69 7377 


Рисунок 8 — Результат уменьшения числа неудачных решений 
с ростом числа циклов обучения 


Для сравнения на рис. 8 горизонтальной линией обозначено количество неудач- 
ных решений для начального варианта системы, полученного на основе экспертных 
оценок. Как видно из приведенного графика, применение СА для задачи адаптации 
системы нечеткого вывода дает лучшие результаты по сравнению с подобранными на 
основе экспертных оценок параметрами на базе априорной информации. Нахождение 
хороших решений происходит за сравнительно небольшое время, что было достигнуто 
путем применения широкого набора оптимизационных техник программирования, в 
том числе и распараллеливания алгоритма. 

Применение генетических операторов. Для формирования новой популяции над 
отобранными хромосомами производится ряд генетических операций. В данной задаче 
применяются два основных оператора — скрещивание (сгоззоуег) и мутация. При этом 
скрещивание происходит при получении 100% новых хромосом. Мутация же произво- 
дится лишь над очень небольшим количеством ген. Оценивание приспособленности 
хромосом в популяции (так же, как и селекция) производится с помощью функции при- 
способленности, определяемой по общему количеству удачных решений для данной 
популяции. 
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они р: 


А.В. Бухановський, Ю.1. Нечаев 

Еволюцйне моделювання задач багатокритер!ально! оптимзацй в систем! 

дослдницького проектування морських суден 

Розглядаеться шджд 1 концептуальний базис еволюшИного моделювання на сучасному еташ розвитку 
вчення про природну генетику в задачах досллдницького проектування. Сформульован! теоретичн! 
принципи формалзаци проблеми з використанням досягнень 1нтелектуальних технолог!й та високо- 
продуктивних засобв обчислень. Наведен! приклади реалзацй розробленого щдходу в практичних засто- 
суваннях дослдницького проектування морських суден. 
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